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Abstract. Based on the nearest-neighbor tight-binding model, we present analyt- 
ical description of low-energy electron transport through a symmetrical junction of 
two semimetal armchair ribbons. The results obtained demonstrate the transmis- 
sion suppression in the vicinity of the neutrality point except for the junction of 
the ribbons consisting of 3p — 1 and 3p + 5 dimer lines, p/2 is odd. Unlike other 
interconnections, these compounds are shown to be free of even local levels arising 
at the junction interface and exhibit electron backscattering to be inversely as the 
square of p + 1. 

Introduction 

Since the discovery of fullerenes, carbon-based materials have been the subject of intense research, 
which led to the discovery of carbon nanotubes and the fabrication of individual one-atom thick graphene 
layers. This opens unprecedented avenues for the investigation of quantum transport in low-dimensional 
ID and 2D systems, as well as attracting the interest of industries, given the potential for innovative 
applications. Motivated by possible device applications, the tr ansport properties o f graphene ribbons 



ap 

(GRs) has been investigated intensive ly both in tight-bindin g IWakabayashi et all 11999c ICresti et al 
20081: iKlvmenko and Sheytsovl 12009a | and Dirac formalisms BeenakkerT 2008 ; Katsnelson et all " 2006 



Blanter and Martini . 120071 ] . Both approaches demonstrate that graphene edges drastically change the 
conducting properties of GRs: zigzag edge produces localized edg e states lead i ng to the metall i c typ e 
of conductivity, while no localized state appears in an armchair GR |Ando . 20051 : Brev and Fertisl . [2006] . 



Within the nearest-neighbor tight-binding model GRs can be either metallic if M — 3p — 1 (p is an 
integer), or semiconducting (M = 3p or M — 3p+ 1). Here M is a ribbon width in terms of dimer lines. 

In this paper, based on the nearest-neighbor tight-binding model, we investigate analytically the 
electron transport through the junctions of armchair graphene ribbons. Specifically, we consider squared 
symmetrical junctions of two semi-infinite graphene ribbons, which have common longitudinal axis (see 
Fig. [2|). Our numerical computations show that there are junctions of a certain type which are almost 
transparent with respect to the incident wave in the vicinity of zero-energy point (K-point). Unlike 
them, the low-energy transmission spectrum of other symmetrical junctions demonstrates suppressed 
transmission. In the present work we find exact analytical expression for the transmission coefficient 
of the junction near the K-point. As follows from the theory, the low-energy behavior of electron 
transmission through a semimetal GRJ is tightly associated with the number C of edge states arising 
at zigzag-shaped boundaries of the ribbons, if they are studied independently. It is important to note, 
that due to the symmetry of the system, transversal modes of different parity in the leads cannot be 
mixed. Thus, since the conducting mode in both channels is even, equal number of even edge states in 
both parts of the junction leads to the practically perfect transmission while a difference between them 
results in the transmission suppression. In the latter ca se, a zero-energy l ocalized state appears at the 



interface between two ribbons. Thus, in agreement with [Li et all l2009bl lal] . the low-energy transmission 
suppression arises from the local level appearance near the junction interface. If a junction is free of the 
even local levels, the transmission exhibits practically perfect propagation. 



Wave solution for electron in armchair GRs 

The infinite armchair GR is considered as a set of elementary cells of two atoms a = A, B with edges 
in armchair and in anti-armchair configurations, see Fig. [TJ We enumerate the cells by two numbers 
{n, m} possessing both integer and half-integer values simultaneously, — oo < n < oo. By taking the 



SHEVTSOV ET AL.: GRAPHENE RIBBON JUNCTIONS 




Figure 1. Sketch of infinite graphene ribbons in armchair configuration (left panel) and in anti-armchair 
configuration (right panel), consisting of M dimer lines in the transverse direction. M is even for anti- 
armchair configuration and odd for armchair one. Black or empty arrows denote cells with integer n, see 
details in text. Shadowed blocks show the nearest neighboring cells. 



tight-binding representation for molecular orbitals we come to the set of linear Schrodinger equations 

-£tyn,ro,A = tpn,m,B + ^ n _i )Tn+ i,B + ^ n _i, m -i,B) ,-js 

with respect to wave function components taken at 2p z orbital of a-th atom in the elementary cell, 
E = -E/I/3I is the electron energy in units of hopping integral \[3\ between the nearest-neighbor carbon 
atoms, (5 < 0. Carbon site-energy is assumed to be zero and serves as the reference. 

Since boundary carbons of infinite armchair GRs have only two neighbors (not three as all inner 
atoms) , it demands the wave function to vanish on the set of absent sites closest to th e boundary 
atoms (Robinson and Schomerusl I2007 ; Zheng et al. , 2007 ; Akhmerov and Beenakkerl . l2008t . Evidently, 
analytical form of the boundary conditions depends on a numeration method of the elementary cells 
along n direction. If index n takes integer value in the regions marked by the black arrows on Fig. [TJ we 
get ipn,o,a = for the lower boundaries, and ipn.N+i.a — 0, ijj n +i/2.N+i/2,a = for the upper boundaries 
in armchair (Fig. [TJi) and anti-armchair (Fig. [TJa) configurations, respectively. Here N = [M/2] is the 
number of elementary cells in the marked regions, M is the number of dimer lines in the transverse 
direction of the armchair ribbon, and [M/2] is the integer part of M/2. Similarly, if we set n to be 
an integer in the regions marked by the empty arrows, one can obtain VVi+i/2,i/2,a = for the lower 
boundaries, and VVi-i/2, N+x/2, a ~ 0, ipn,N+i,a = for the upper ribbon boundaries on the corresponding 
Figures. Here N = [M/2] + 1 for armchair configuration and N = [M/2] for anti-armchair one. 

The nontrivial solutions to Eq. JT|) satisfying the appropriate boundary conditions can be written as 
follows 

,/. _ XjM „ifc„„ f 1, a = A, 



VM+T 



a = B. 



(2) 



Here we introduced the transversal wave functions for the numeration methods marked by black and 
empty arrows respectively 



Xj (m) = sin(^m), Xj{ m ) = sin fe( m - 1/2)], 



(3) 



where quantity 



27TJ 

m + r 



i<i< 



M 

T 



(4) 



has meaning of the dimensionless transversal wave number for electro n in armchair GRs. The longitudinal 
wave number kj from Eq. ([2]) can be found from dispersion relation |Andol . 2005 1 . which can be derived 
by plugging in Eq. @ to Eq. (fTJ) and demanding zero determinant. Thus, 
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and quantity 9j is denned by expression 



l + 2coslf-e i ^ /2 
E ' 



(0) 



Given E and £j, real- valued number kj determines a propagating state in the channel. Here 
real-valued quantity 0j has mea ning of phase shift between A and B atoms in the elementary cell 
Klvmenko and Shevtsov , 2009a and state ([2]) is normalized on a unit celfj] of the ribbon. The group 



velocity connected with the propagating state can be written as follows 



dE 



1 £j . kj 
- — cos -f sin ~~~ 



E 



Thus, the propagating wave ((2]) runs from left to right in the valence subbands and vice versa for the 
conducting ones. 

It should be noted that the choice of the elementary cell numeration along n direction of the infinite 
graphene ribbons remains arbitrary. Obviously, the value of the wave function component ip n ,m,a does 
not depend on certain choice. Considering armchair GR junctions, we'll refer elementary cells n — to 
the junction interface. In this case, the numeration of elementary cells in the ribbons becomes unique. 

General equations for scattering amplitudes 

Now we consider a symmetrical junction of two semi-infinite graphene ribbons with armchair edges, 
which have common longitudinal axis (see Fig. [2|). Such junctions can be, obviously, constructed only 
of ribbons in armchair configuration. The corresponding matching conditions are as follows 



0, 1 < m < /i, 

%, m ,A = { ^o, m .A » H+l<m<Ni + n, 
0, Nl 1 < to < N r . 



(7) 
(8) 



Later on, by indexing c = I,r we will distinguish quantities referring to the corresponding left (I) and 
right (r) leads. 



m=Nl+\l 




right 



left right 



Figure 2. Two possible geometries of the symmetrical junction of armchair GRs. Shadowed rectangle 
depicts the scattering region. Transition from case a to case b is provided by deletion of one zigzag line 
from the shoulders, and vice versa. Figure represents junctions with three dimer lines long shoulders as 
they are the subject of the main text. 

The incident wave with the transverse mode jo in the left lead, meeting the interface between 
two ribbons, scatters into both propagating and evanescent modes. The general wave solutions for the 



1 By definition, the unit cell of the ribbon is a union of elementary cells along the translation period of the ribbon. It 
consists of all {n, m} and {n + 1/2, m} elementary cells, where the number n is fixed. 
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electron in both leads should be constructed as a linear combination of expression ^ over all transverse 
modes j in the corresponding channel. In particular, for E < we obtain 



[Mi/2] [M r /2] 

n,m,a 



(9) 



Here co-factors rj j0 and fjj for the propagating mode j have meaning of the scattering amplitudes in 
the corresponding channel. When fcj is complex, the state ipn,m,a(j) describes an evanescent eigenmode 
in the corresponding lead. In this case r 3 - J0 and tjj are unphysical quantities and there is no need to 
normalize C,m,a(i)- 

Plugging in expressions ^ into conditions ([7]) , ijSJ) and taking into account orthogonality of functions 
{x^(to)}, we come to a set of equations 



[Mr/2] 

x i -iO l - \ " i8 r ,, 

e 3 °Mo+ e 3r J,jo= 2^ 9j,j' e 3t f,jo 
j'=l 

[Mi [2] 

tj,jo ~ 9j',j [Gj',jo ~^~ r j',jo\ ' 

i'=i 



(10) 



where 5^ J0 is the Kronecker delta and 



9j,f 



Ni 



y/(M t + 1)^ + 1)^ 



= Z 4(m)xy(m + /i) 



(11) 



is a matrix coupling transverse modes of the left and right leads. 

Eqs. (|10p hold for arbitrary types of ribbons (scmimctal or semiconducting) forming a squared 
junction with a common edge. Solving them with respect to scattering amplitudes tj j , one can obtain the 
transmission probabilities Tjj — (vj /Vj )\tj^ \ 2 and the transmission spectrum T(E) = ■ ■ Tjj [E) 
of the junction. 




-case a, Mr 17, M,=23 
-case b, Mr 17, M,=23 



-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 




Figure 3. Transmission spectrum for semimetal 17/23 and 23/29 graphene ribbon junctions. The 
corresponding junction geometry is specified on Fig. [2 where cases a and b refer to the left and right 
panels correspondingly. 

Later on, we focus on the metallic graphene junction Mi/M r where M v = 2>p v — 1, and p v is an 
integer. Numerical modeling for 17/23 junctions, the left panel on Fig. [31 demonstrates the nonzero 
transmission, however, calculations carried out for 23/29 junctions show transmission suppression in the 
vicinity of zero-energy point. Thus, we can find that junction family (3p — l)/(3p + 5) with odd p/2 
demonstrates practically perfect transmission while the transmission is suppressed for even p/2. More 
detailed analysis of transmission properties of semimetal graphene junctions observed at E — is the 
subject of the next section. 

Propagation at the neutrality point in semimetal junctions 

Detailed analysis of band structure based on Eq. J^J shows that in the vicinity of zero-energy point 
(E = 0), the only transversal mode available for propagating has number j v — p v in each channel. 
Meantime all other modes are vanishing since in that energy region Eq. Jl| does not give real-valued 
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longitudinal component of the wave vector, kj. By substitution of £^ = 2tt/3 to energy dispersion ([5]) 
we obtain k pi/ = 2ir. To find 9 pv we plug in the energy dispersion to expression for cos 6" following from 
Eq. ([6]), 9 pi/ = 7r/2. Thus, one can write down the following expression for the propagating waves in the 
leads 

VnmM= if** , i •' ""r V = l,r. (12) 

Evanescent electron state in the leads corresponds to a complex- valued fcj. It can also be constructed 
from state (J2j) by substituting = 7r + iSj /2 into definition (O taken for E = 0, 



2 COS -y , 



2cos^- > 1 <^ j <p„, 



l/[2cos^], 2cos% < 1 < j < [^] 

With help of Eqs. (O and (fT5]) we must choose regular evanescent solutions for both leads. Finally. 



(13) 



1pLm,cx(J >Pl) 



v / MTFT 
_ X'jjrn) gl 



Yn,m,a\J ^FIJ + 1 




Yn,m,a\J ^ Pr ) — ^ M r + 1 



0, 


a 


-A, 


1, 


a 


= B, 


1. 


a 


= A, 


0, 


Q 


= B. 



(14) 



~ M v ~ 




Pu — 1 


2 


~Pu = 


2 



It is important to note that the evanescent eigenmodes j > p v in Eqs. (|14[) coincide identically with 
the edge state solutions arising at zigzag-shaped boundaries of the corresponding isolated semi-infinite 
leads wh ere the appropriate boundary condition s ij) l n n (j > pi) = 0, V'StoaO > Pr) — must be 
satisfied Wa kabavashi et al. . 19991 : Nakada et al. , 19961 ] . The number of such modes (with j > p v ) in 
the corresponding ribbon is determined by dependence 



since 1 < j < [^] . 

To obtain equations for scattering amplitudes, one needs to satisfy interface conditions ([7]) and l[8|), 
using the linear combination of solutions (|12|l and (fT"4|) . The result reads, 

Pr + Cr 

i5 ],p l ( 2 ~ f J.,Pi)+ f 3,Pi r l(Pl -3)= 9j,j'tj>ml i§ 3',Pr +V(f ~Pr)l 

tj,Pl Ihvr + ViPr - j)] = 9j'.A 5 3'-,Pi r J,Pl + r f,PlV(j' - Pi)]- 

Here f J;Pl = rj tPl + Sj iPl , r)(x) = 1 if x > and rj{x) — otherwise. Using Eq. ([3]) and definition (H|), one 
can find the following properties of the coupling matrix gjj/ . First, due to the symmetry of the system, 
transversal modes of different parity cannot be mixed, which means that gjji =/= only if j and j' are 
odd or even simultaneously. Second, when j' = p r , one can verify that 



9j,p r = 9$j,Pi , 9 = y/pi/ipi + 2 )- 



(16) 



Since the propagating mode in both leads is even, one must solve Eqs. (|15[) taking into account only 
even-numbered modes. This is the subject of the next section. 



Analytical solution for transmission probability 

First we introduce the number of even edge states in the leads, £° vcn = [£„/2]. Due to relation P^]) . 
the second equation in (fT31) taken for j > p r is simplified to the set of linearly independent equations 



9pi ,j r Pl ,Pl + 

j> -p, =2,4, 



tprtPl 9 T Pl ,Pl 
9j',3 r 3',Pl = 0, j — Pr = 2,4, 



/*cv 



(17) 
(18) 



If £E 



> Cf 



- 1, homogeneous system (fT5| always has trivial solution Tj> Pl , Pl = 0, which leads to the 
transmission suppression T = \t Pr . Pl | = 0. If C° r vcn = £f vcn , then the first equation in (fT5l) taken for even 
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j > pi leads to the closed set of homogeneous equations. That system has trivial solutions tj >Pr ^ Pl = 0. 
Substitution of this result to the first equation in (1151) taken for j = pi gives the equation 



9^p r ,pi ~r" r pi,pi ~~ 2. (19) 

Solving Eqs. (fl~7f and (|19p with respect to the transmission amplitude t Pr , Pl , one can obtain the following 
expression for transmission coefficient 



Condition £° ven = £° von holds only when p r = pi + 2 and pi/2 is odd. Thus, we obtain 

T=1-( P; + 1)- 2 . (21) 

There is a physical reason for such behavior of the transmission coefficient. Forming a junction be- 
tween two ribbons may lead to local level occurrence near the junction interface. To obtain the condition 
for existence of the states localized at the junction interface one should use only linear combination of 
evanescent solutions (fT4|) omitting the propagating ones. Utilizing boundary conditions ([7]) and ([8]) one 
can find the linear homogeneous system of equations. It has nontrivial solution only when Cf cn ^ £° ven . 
Thus, transparent junctions mentioned above are free of even local levels which causes good conducting 
properties of such systems. In contrast to them other junctions have zero-energy local level of degeneracy 
V = C c r vcn — £° vcn which leads to the full suppression of transmission in the vicinity of E = 0. 



Conclusions 

The results obtained in the present work show that the low-energy propagation through a sym- 
metrical graphene junction is directly connected to existence of even localized states near the junction 
interface. If junctions are free of such local levels, they exhibit almost perfect propagation described 
by dependence (|2"T|) . Other GRJs possess local level suppressing the propagation in the vicinity of the 
K-point. 

Symmetrical junctions do not exhaust all types of transparent junctions. There are also common- 
edged (3p— 1) / (3p+2) (p is odd) junctions with three dimer lines long shoulder which are free of local levels 
too. Applying the same technique as for symmetrical junctions one can find tha t for common-edged junc- 
tions t he reflection coefficient turns out to be inversely as the square of (2pi + l) [Klvmenko and Shevtsov . 
l20Q9b| . Our calculations performed with help of Eq. (flT)]) and Eq. (TIT]) for such systems (specifically, 
20/23 and 23/26 junctions) are i n perfect agreem ent with those obtained numerically in the framework 
of Landauer-Buttiker formalism [Li et all l2009bl la| . The calculations made show that common-edged 
junctions with upward-shifting the narrower lead on one or two dimer lines conserves the unsuppressed 
transmission, though the solution for the transmission coefficient cannot be found analytically. 

Obtained results give an opportunity to go further and consider more complicated systems, such 
as: quantum dot structure, multiple quantum dot structures, superlattices, constrictions of armchair 
graphene ribbons, junction of armchair graphene ribbon and zigzag nanotube, etc. 
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